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Licensing provisions: none 

Computer for which the program is designed and others on which it has been 
tested: 

Computers: (i) IBM PC, (ii) cluster 

Installations: (i) Brock University, St. Catharines, Ontario, Canada L2S 3A1; 
(ii) SHARCNET http://www.sharcnet.ca 

Operating system under which the program has been tested: LINUX 

Programming language used: REDUCE 3.7, REDUCE 3.8 

Memory required to execute with typical data: problem dependent (10 Mb - 
1Gb), typical working size < 100 Mb 

No. of bits in a word: 32, 64 

No. of processors used: (i) 1; (ii) multiple 

Peripherals used: no 

Has the code been vectorizedl no 

No. of lines in distributed program including documentation file, test data, 
etc.: 2485 

Distribution format: ASCII 

Online access and tutorials: http://lie.math.brocku.ca/crack/susy/ 

Keywords: Integrable super-systems; symmetries; recursions; classification; sym- 
bofic computation; REDUCE 

Nature of physical problem: The program allows the classification of iV > 
1 supersymmetric nonlinear scaling-invariant evolution equations {ft = (fi^ , 
bt — (fi''} that admit infinitely many local symmetries propagated by recursion 
operators; here b{x, t; 9) is the set of bosonic super-fields and f{x, t; 9) are 
fermionic super-fields. 

Method of solution: First, (half-)integer weights |/|, . . ., I^tl, \Dx\ = 1 are 
assigned to all variables and derivatives and then pairs of commuting fiows 
that are homogeneous w.r.t. these weights are constructed. Secondly, the seeds 
of higher symmetry sequences [2] for the systems are sorted out, and finally 
the recursion operators that generate the symmetries are obtained [3]. The 
intermediate algebraic systems upon the undetermined coefficients are solved 
by using [4]. 
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Restrictions on the complexity of the problem: Computation of symmetries of 
high differential order for very large evolutionary systems may cause memory 
restrictions. Additional size/time restrictions may occur if the homogeneity 
weights of some super-fields are non-positive, see section 1.2 of the Long Write- 
Up. 

Typical running time: depends on the size and complexity of the input system 
and varies between seconds and minutes. 

Unusual features of the program: SsToOLS has been extensively tested using 
hundreds of PDE systems within three years on UNIX-based PC-machines. 
SsToOLS is applicable to the computation of symmetries, conservation laws, 
and Hamiltonian structures for A'^ > 1 evolutionary super-systems with any N. 
SsToOLS is also useful for performing extensive arithmetic of general nature 
including differentiations of super-field expressions. 
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LONG WRITE-UP 



Introduction 

The principle of symmetry belongs to the foundations of modern mathematical 
physics. The differential equations that constitute integrable [1] models practi- 
cally always admit symmetry transformations and, reciprocally, new classes of 
integrable phenomena in physics are obtained by postulating some symmetry 
invariance. The presence of symmetry transformations in a system yields two 
types of explicit solutions: those which are invariant under a transformation 
(sub)group and, secondly, the solutions obtained by propagating a known so- 
lution by the same group. The two schemes for constructing new solutions of 
PDE are crucial for systems of super-equations of mathematical physics (e.g., 
supergravity models); these equations involve commuting (bosonic, or 'even') 
and anticommuting (fermionic, or 'odd') independent variables and/or un- 
known functions. Indeed, other methods for solving nonlinear equations need 
a special adaptation to the super- field setting, see e.g. [2]; another approach to 
integrability of supersymmetric equations has been investigated in [3] . Hence 
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the symmetry considerations [1], which arc based on the compiitation of in- 
finitesimal symmetry generators for PDE, become highly important. The aris- 
ing computational problems are unmanageable without computer algebra that 
permits handling relevant systems. 

In the literature it has been observed (see [1,4] and references therein) that 
the principal phenomena in nature arc governed by systems of PDE that 
admit higher symmetries, that is, the symmetry transformations that involve 
higher order derivatives of the unknown functions. Simultaneously, when deal- 
ing with supersymmetric models of theoretical physics, it is often hard to pre- 
dict whether a certain mathematical approximation will be truly integrable or 
not. Therefore we apply the symbolic computational approach to the physical 
problem of classifying the systems that exhibit necessary integrability fea- 
tures. It must be noted that the very idea to filter out integrable cases using 
the presence of 'many' symmetries is widely accepted in the computer branch 
of modern mathematical physics, see e.g. [4,5]. These systems are called sym- 
metry integrable [1,4]); in some cases, they can be transformed to exactly 
solvable equations or their extensions. In view of this classification task, we 
analyzed fermionic extensions of the Burgers and Boussinesq equations and 
related the former with evolutionary systems on associative algebras in [6] . 

This paper is organized as follows. First we formulate the axioms of the clas- 
sification problem for = 1 supersymmetric systems of evolutionary PDE. 
In sections 1.1 and 1.2 we describe the two modes of the procedure ssym 
in the package SsTooLS that allow, respectively, finding unordered pairs of 
commuting fiows and the computation of symmetries for previously found sys- 
tems. Also, in section 1.2 we review a geometric (coordinate-free) algorithm 
for constructing recursion operators. However, we refer to [5,7] for basic no- 
tions and concepts in the geometry of (super)PDE, see also [8,9,10,16] and 
references therein. The principal result of this paper is that there exist only 
four nonlinear coupled boson-fermion systems that satisfy the axioms and the 
weighting |/| = |6| = \Dt\ — \. Section 1 concludes with the classification of 
their recursion operators. In the following sections we investigate properties 
of these four systems, giving simultaneously the syntax and describing appli- 
cability of SsToOLS subroutines developed for the calculation of the scaling 
weights, symmetries, linearizations, conservation laws needed for introducing 
the nonlocalities, and recursion operators. Finally, sample runs of SsToOLS 
are given. 



1 The classification problem 

Let us introduce some notation. We denote by 6 the super-variable and we put 
V = De + OD^ such that = and [V,V] = 2D^; here Dq and are the 
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derivatives w.r.t. 9 and x, respectively. Throughout this text, the operator V 
acts on the succeeding super-field unless stated otherwise explicitly. 



The package SsToOLS in Reduce (ver. 3.7 or higher) is designed for calcu- 
lus on super-PDE. It was implemented for solving the classification problem 
for symmetry integrable > 1 super-PDE and fermionic extensions of the 
purely bosonic evolution equations. Compared to other programs listed in the 
overview [11], it is able to work with (super)PDE and compared to programs 
mentioned in [12,13,14] the package SsToOLS has not only the ability to com- 
pute higher symmetries for a given system but also to solve the non-linear 
problem of computing polynomial systems and their higher symmetries at 
once if certain homogeneity weights are specified; here we also note the pro- 
gram [15] for computations performed in presence of the supersymmetry. In 
SsToOLS, the following list of seven axioms was postulated: 



(2) all equations are spatial translation invariant and do not depend on the 
time t explicitly; 

(3) none of the evolution equations involves only one field and hence none of 
the right-hand sides vanishes; 

(4) at least one of the right-hand sides in either the evolution equation or its 
symmetry is nonlinear; 

(5) at least one equation in a system or at least one component of its sym- 
metry contains a fermionic field or the super-derivative V; 

(6) the equations are scaling invariant: their right-hand sides are differential 
polynomials homogeneous with respect to a set of (half-) integer weig htB 
\V\ = i, \Dx\ = 1, \Dt\ > 0, |/|, \b\ > 0; we also assume that the positive 

weight of Ds that occurs through (J^"^ is (half-) integer. 

The time t and the parameter s along the integral trajectories of the symmetry 
fields can either assign the bosonic fiows bt, hg to h and the fermionic fiows ft-, 
fs to / or reverse the parities, that is, bt or 6j become fermionic and ft or /j 
are bosonic. Here we use the notation t or s if the parities of ft, bt or /§, bg 
are opposite to the parities of / and b, respectively. 



^ The weights \ ■ \ take all variables (e.g., a fermionic super-field /) and deriva- 
tions (e.g., Dt = d/dt) to (half-)integers (|/| and \Dt\ € ^Z, respectively) such 
that the weight of any product is the sum of weights of the factors (hence 
\ft\ = \Dt{f)\ = \Dt\ + I /I); only homogeneous algebraic expressions with equal 
weights of the components are considered from now on. 




admits at least one higher symmetry 




5 



1.1 Computation of commuting flows 



The procedure ssym for computing symmetries, which is the central proce- 
durqfj of the SsToOLS package, can be used in two modes. The first mode 
of ssym developed together with W. Neun was used for finding N = 1 super- 
symmetric coupled boson-fermion evolutionary systems that satisfy the above 
axioms. The axioms are incorporated in the ansatz generators. Hence, having 
the weights |/|, \t\ or and |s| or |s| specified, we face the problem of 
constructing a pair of commuting flows (/, &)t(or t) and (/, &)s(or s)- The proce- 
dure ssym generates the homogeneous ansatz with undetermined coefficients 
in both flows, and calls the program Crack [16] for solving automatically the 
overdetermined algebraic systems. The standard call in the first mode is 



ssym(N , tw , sw , af wlist , abwlist , eqnlist , f 1 , inelist , flags) ; 



where 




N 


. . the number of superfields 


tw 


. . 2 times the weight of Dt\ 


sw 


. . 2 times the weight of Dg] 


afwlist . 


. . list of weights of the fermion fields f (1), f (2),. . ., f (nf ); 


abwlist . 


. . hst of weights of the boson fields b(l), b(2),. . ., b(nb); 


eqnlist . 


. . list of extra conditions on the undetermined coefficients; 


fl 


. . extra unknowns in eqnlist to be determined; 


inelist . 


. . a list, each element of it is a non-zero expression or a list with 




at least one of its elements being non-zero; 


flags 


. . init: only initialization of global data. 




zerocoeff : all coefficients = which do not appear in inelist 




tpar: if the time variable t changes parity. 




spar: if the symmetry variable s changes parity. 



The nonnegative integer N is the number of super-fields 9^ such that Vi = 
Dgi + 9'' Dx are the respective super-derivatives, Vf = D^. The computer 
notation for is d(i , . . .) , thus opening the opportunity to work with > 1 
under SsToOLS doing the calculus of super-fields; here we also note that 
the syntax of is df (. . ,x) and similarly for Dt. In addition, SsToOLS 
permits the investigation of non-supersymmetric systems (with = 0, which 
is not needed for the four systems within the classification below). What may 
be of interest to the specialists working with super-calculus and super-PDE 
is that after loading SsToOLS the simplification of polynomial expressions 

^ The command sshelpO ; provides a detailed user's guide on the package includ- 
ing examples that illustrate admissible combinations of the flags in the procedure 
calls. All package subroutines are described in this paper where appropriate. On- 
line access to SsToOLS without any need to install the computer algebra system is 
provided through http://lie.math.brocku.ca/crack/susy. 
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which involve the anti-commuting derivations Vi and fields p is performed 
automatically, taking care of all the minus signs that arise in multiplications 
and derivatives. 

The numbers nf of fermion fields and nb of boson fields are determined auto- 
matically through the number of elements in the lists afwlist and abwlist. 
The input list eqnlist can contain substitutions, like p2=l, or expressions 
which are to be set to zero, like p3*r4+p2*r3. Conditions as p2=l are exe- 
cuted instantly at the time of formulating the ansatz for the pair of commuting 
fiows, while expressions without = sign are added as conditions . . .=0 to the 
other equations when calling Crack. Typically, one would run the program 
first with the flag init to see which ansatz for the system and its symmetry is 
generated and then start ssym again without init but with the option to add 
extra conditions on the unknown coefficients either through the flag zerocoef f 
or through specific extra conditions in eqnlist or entries in inelist. 

Using the bounds < |/|, |6| < 5 and < \Dt\ < \Ds\ < 5, we constructed the 
experimental database [17], which contains 1830 equations such that N < 2 
(the duplicate PDEs that appeared owing to possible non-uniqueness of the 
weights are not counted) and their 4153 symmetries (plus the translations 
along X and t, and plus the scalings whose number is in fact infinite). In this 
paper we investigate properties of the systems for which the weights of the 
times t, i are \t\ = \t\ — — | 4^ \Dt\ = \Dt\ = |, assuming further that 
I/I ~ 1^1 ~ i (^^® weights may not be uniquely defined). Having fixed the 
weights and using SsToOLS, we discovered four systems that satisfy the set 
of axioms; they are classified as follows, see Table 1 on p. 12 below. System (2) 
is related to the fermionic extension of the Burgers equation [6]. Also, we get 
a super-field representation (4) of the Burgers equation (5) and discuss its 
properties which are revealed only by the presence of super- variables. Further, 
there is a multiplet of 14 systems (9) with the parity preserving time t; these 
systems possess s-symmctries and higher s-supersymmetry flows. Finally, we 
obtain a unique system (11) which has the parity reversing times t and s; using 
this system, we illustrate the process of calculating the recursion operators in 
the test output. 



1.2 Reconstruction of recursion operators 

The second application of SsToOLS is the construction of symmetries, con- 
servation laws, and recursion operators for the symmetry algebras of given 
systems, i.e. systems S that have been obtained before by ssym and placed in 
the database [17]. The method of Cartan's forms [5] for the recursion opera- 
tors is applied, see [6] for examples based on the notation introduced below. 
Within this approach, the recursions TZ are regarded as symmetries = TZ^ , 
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Bgj^ = TZ'^ of the linearized equations Lin£^. Namely, we ignore the differential 
function structure of the symmetry flows = F , hg = B and regard F and 
B as the components of solutions of the linearized equations Lin£, which are 
defined in the next paragraph. The expressions TZ — R{F, B) are recursion 
operators for £ if each TZ satisfies the linearized equation Lin£^ again and if 
they are linear w.r.t. F, B, and their derivatives. 

The understanding of linearized systems Lin £ from a computational viewpoint 
is as follows; we consider the differential polynomial case since this is what 
we and SsToOLS deal with. Given a system formally assign the new 'lin- 
earized' fields F' = f (nf+i) and B^ = b(nb+j) to = f (i) and 6^' = b(j), 
respectively, with i = 1,. . .,nf and j = 1,. . .,nb. Pass through all equations, 
and whenever a power of a derivative of a variable /* or b' is met, differen- 
tiate (in the usual sense) this power with respect to its base, multiply the 
result from the right by the same order derivative of or B^ , respectively, 
and insert the product in the position where the power of the derivative was 
met. Now proceed by the Leibnitz rule. The final result, when all equations in 
the system £ are processed, is the linearized system Lin£^. For example, the 
linearized counterpart of bt = b{Vff is Bt ^ B ■ {Vff + 2hVf ■ VF. 

The linearization Lin £ for a system of evolution equations £ is obtained using 
the procedure linearize: 

linearize (pdes, nf ,nb) ; 

where 

pdes . . . list of equations with first order derivatives in l.h.s.; 
nf ... number of the fermion fields f (1) , . . .,f (nf ) ; 
nb ... number of the boson fields b ( 1) , . . .,b (nb) . 

The values nf and nb may be greater than one for the coupled boson-fermion 
systems with unique / and b (only such systems are studied throughout this 
paper) because some additional variables (the nonlocalities, see (1) for a de- 
tailed explanation) can have been added previously by hand using the method 
of trivializing conservation laws. The actual number of fields is doubled by a 
call of linearize, and in the sequel we always denote these new 'linearization 
variables' by the respective capital letters. 

Now we describe the syntax of a call in the second mode of using procedure 
ssym. In this mode symmetries are computed (which can be interpreted later 
as recursion operators) for a given evolutionary system pdes. The call is 

depend {f (1) ,f (2) , . . ,f (nf ) ,b(l) ,b(2) , . . ,b(nb)},x,t; 
max_deg : = . . ; 

hom_wei:={{sw_l,afwlist_l,abwlist_l}-, . .}; 
ssym(N,tw, sw, af wlist , abwlist ,pdes ,f 1 , inelist , flags) ; 
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where this time we use the actual numbers of fields (they may have increased 
after linearize produced a bigger system by assigning new to /* and 
to 6^). We now have 
pdes . . . the list of equations df (. . . ,t)=. . . specifying the (linearized) 
system for which symmetries are to be computed, as well as 
substitution rules . . .=>. . . and other conditions ...=... (see Re- 
mark 1 below); 
flags . . . the list of flags now includes: 

lin: the symmetry to be determined is linear w.r.t. the fields 

f (i), b(j) such that i > nf/2 and j > nb/2; 

filter: if present this flag indicates that the homogeneities 

listed in hom_wei have to be used in addition; 

tpar, spar: and Dg, respectively, is parity changing. 

If the right-hand sides in pdes involve any constants which should be computed 
such that the higher symmetries do exist, then these constants are listed in the 
input list fl, like {pl,p2,. . .} otherwise fl is {}. This may be needed when 
we suspect the presence of a multiplet of systems parameterized by discrete 
admissible values of these coefficients pl,p2,. ... If the constant coefficients are 
not given in the list f 1, then they are treated as independent parameters and 
only symmetries (and not these constants) are determined for generic values 
of the constant coefficients. 

Next, if a boson's weight is non-positive or a fermion's weight is negative 
then the global variable ina:x_deg must have a positive integer value, which is 
the highest power of such fields and their derivatives in any ansatz that the 
program will generate. 

Remark 1 The key for procedure ssym to know whether to compute sym- 
metries for the given system (this is mode 2) or to make an ansatz for the 
yet unknown system with a symmetry (hence, for a pair of commuting fiows 
within mode 1) is the first element of the input list pdes. If this is of type 
...=... or of type . . .=>. . . and if the left hand side is df (. . .) or d(. . .) , then: 

• All elements of pdes of the form df (. . . ,t)=. . . are collected and interpreted 
as the system for which symmetries are to be computed. 

• AH other elements of pdes of the form . . .=>. . . are assumed to have df (. . .) 
or d(...) on the left-hand side. These substitution rules and differential 
consequences of them are used to simplify right-hand sides of the system, 
in making the ansatz for the symmetry, and in further computations. 

• All other elements of type ...=... or elements that are simple expressions 
(which are set to 0) are interpreted as extra conditions on the undetermined 
coefficients. 

There are two types of relations to be provided in pdes which is the reason 
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why substitutions => are needed in addition to the usual relations =. First, the 
very algorithm [5,10] for finding recursions of symmetry algebras of equations E 
suggests to treat the recursions as symmetries of the linearized systems Lin£, 
taking into account the original equations (that is, replacing their left-hand 
sides with their right-hand sides) and differential consequences of them. The 
linearized system, which involves time derivatives of Cartan's forms, is in- 
cluded in the list pdes in form of relations df (f (i) ,t)=. . ., % > nf/2 and 
df (b(j),t)=. .., J > nb/2. 

Secondly, we recall that the recursion operators are frequently nonlocal, thus 
requiring introduction of new nonlocal variables for the systems S, see [5,6] 
for many examples. To this end, the systems are supplemented by differential 
relations (see (1) in section 1.3) which specify the new nonlocal dependent 
variables. Then these additional relations and their linearized counterparts 
are substituted but their symmetry flows are discarded, which is equivalent 
to setting them to zero. These substitution rules are included in pdes us- 
ing => in the form df (f (i) ,x)=>. . ., d(l,f (i))=>. . . and df (b(j) ,x)=>. . ., 
d(l,b(j))=>.... 

The order in which the = and => relations occur in pdes is inessential. 

Suppose the system in pdes is homogeneous w.r.t. several sets of homogeneity 
weights simultaneously, that is, it admits several scaling symmetries at once. 
One then wants to find a symmetry that obeys the same sets of homogeneity 
weights for the super-fields, each set supplemented by the expected weight of 
the parameter s. This is done as follows. An arbitrarily chosen (see a comment 
below) set of the weights tw,sw for the times t,s and af wlist , abwlist for 
the super-fields fi, bj is specified among the arguments of ssym, while the 
other weights are given through the global variable hom_wei. This is a list of 
lists {sw_i,afwlist_i,abwlist_i} with the additional sets of homogeneity 
weights of the parameter s and the dependent variables, respectively. The fiag 
filter causes computation only of those symmetries that satisfy all additional 
homogeneities encoded in the global variable hom_wei as well. 

The choice of the primary set of weights placed in the ssym arguments is ar- 
bitrary with the following comment: if all bosons' weights arc positive and all 
fermions' weights are nonncgative in the primary set, then the bound max_deg 
is not applied even if hom_wei contains zero or negative weights of any vari- 
ables. Indeed, in this case the primary weights already guarantee a finite num- 
ber of terms in the ansatz for the symmetry, and the extra homogeneities in 
hom_wei act as an additional filter. The global bound max_deg is needed if for 
each set of homogeneity weights either one boson's weight is non-positive or 
one fermion's weight is negative. 

Remcirk 2 The (non) uniqueness of the weights for a given system can be 
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checked within some range by using the procedure wgts: 

wgt s (pde s , nf , nb , maxwt ) ; 
where 

pdes . . . hst of equations df (. . . ,t) = . . .; 

nf ... number of the fermion fields f (1),. . .,f (nf); 

nb ... number of the boson fields b(l),. . .,b(nb); 

maxwt . . . only weights with a total sum < maxwt are to be listed. 

The weights are scaled such that the weight of D^. is always 1 and the weight 
of P is |. Note that the computer representation of the weights is twice the no- 
tation used throughout the paper to avoid half-integer values. This convention 
is usual. 

Remark 3 The systems that admit several scaling symmetries and hence are 
homogeneous w.r.t. different weights allow to apply the breadth search method 
for recursions, which is the following. Let a recursion of weight \sr\ w.r.t. a 
particular set of weights for the super-fields /, h and the time t be known. 
Now, recalculate its weight |s^| w.r.t. another set of homogeneity weights |/|, 
|6| and then find all recursion operators of weight \s'^\. The list of solutions 
will incorporate the known recursion and, possibly, other operators. Generally, 
their weights will be different from the weight of the original recursion w.r.t. 
the initial set of weights. Hence we repeat the reasonings for each new operator 
and thus select the weights \sr\ such that the recursions exist. This method is 
a serious instrument for checking consistency of calculations and elimination 
of errors since the existence of a solution is always guaranteed. We used it 
while testing the second mode of ssym in the SsToOLS package. An example 
of implementing it is given in the test output, where the recursion (10) for the 
multiple homogeneous a — 2 system (9) is constructed. 



1.3 The classification of systems and recursions 

Let us introduce some notation which will simplify classification of recursion 
operators for the systems at hand. Assume TZ is a. recursion for an equation 
and consider the symbol OTcT'^'^'^leight • The subscripts 'ord' and 'weight' denote 
the differential order and the weight of the recursion 7^, respectively. 

Recall further that the nonlocal variables are defined for = 1 super-equati- 
ons £ by trivializing [5,6,9] conserved currents, which are of the form Dt{p) + 
VQ = 0; here the right-hand side vanishes by virtue (=) of the system S 
and all possible differential consequences from it. The standard procedure [5] 
suggests that this relation is the compatibility condition for a new 'nonlocal' 
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variable (the nonlocality ) , say w, whose derivatives are set to 

Vw := p and Wt := ^Q, (1) 

where the sign corresponds to the parity preserving ('— ') or reversing ('+') 
time t and t, respectively. Each nonlocality thus makes the conserved current 
trivial; the new variables can be bosonic or fermionic. Hence, starting with 
an equation one calculates several conserved currents for it and trivializes 
them by introducing a layer of nonlocalities whose derivatives are still local 
differential functions. This way the number of fields is increased and the system 
is extended by new substitution rules. Moreover, it may acquire new conserved 
currents that depend on the nonlocalities and thus specify the second layer 
of nonlocal variables with nonlocal derivatives. Clearly, the procedure is self- 
reproducing. An example of fixing the derivatives of a nonlocal variable is given 
in the text output. So, one keeps computing conserved currents and adding 
the layers of nonlinearities until an extended system £ is achieved such that 
its linearization Lin£^ has a (shadow of a, [5]) symmetry TZ\ this symmetry of 
Lin£ is precisely the resursion for the extended system £. 

The superscript 'layers' (if non-empty) in the symbol OTcT'^^^^weight indicates 
the required number of layers of the nonlocal variables assigned to conserved 
currents. The symbol 'y' denotes the number of recursions found for a given 
differential order, weight, and given nonlocalities. In the sequel, we denote 
local recursion operators by L and nonlocal ones by N. The symbol Z is used 
to denote a nilpotent recursion operator whose powers are equal zero except 
for a finite set of them, and is a super-recursion that swaps the parities of 
the flows. 

Tabic 1 

The classification of coupled super-systems and their recursions with respect to the 



primary weights |/| = \h 


= lAl = lA-l = i 


(2) 1 


ft = Vb + fb, 
ht = Vf 


lN\„\N^_„lN\,,lNi, 

2 2 2 2 2 


(4) 1 


ft = -Db, 

bt = Vf + b^ 


l/\rl 


(9) 1 


ft = -afb, 
_ bt = Vf + b^ 


0. rl rl O7I Oynl Oynl Oyl 
a — Z. l-£^_25 l-L' nl, 0-^-2' 0^-25 0^ 9I' 0^ 9I' 
2 2-^2 2 2 


(11) 


{fj = Vf + b\ 
\bj = Vb + fb 


qU^ 1 , iL^i, 3 
-52 2-2 



It turns out that the equations presented in Table 1 exhibit practically the 
whole variety of properties that super-PDE of mathematical physics possess. 
Let us discuss them in more detail. 
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2 A fermionic extension of the Burgers equation 

The system 

ft^Vb + fb, bt^Vf (2) 

is homogeneous w.r.t. a unique set of weights |/| = \b\ = |, \Dt\ = |, = 1. 
System (2) admits symmetries (/s, bg) for all weights \Ds\ > |. 

We recognize that system (2) is related to the fermionic extension of the 
Burgers equation [6] . Namely, consider the fermionic super-field = r]{t, x) + 
6$,{t,x) of weight whose derivatives are 0t = /, T>(f) = b and which thus 
potentiates the second equation in (2). Then from the first equation in (2) 
we get the scalar equation (p^t = 4>x + (pt ' T^'P- Its component form is the 
evolutionary system 

= mt - im, ix = itt -iit + oi- ritrju, a = 1, (3) 

where the variable x is the evolution parameter (the time) and the coordinate t 
is the new spatial variable. System (3) is precisely the a — 1 fermionic (su- 

per)extcnsion of the Burgers equation upon the functions b{t,x) and w(t,x), 
see [6, Eq. (16)], rewritten by setting ^{t,x) = —b{x,t) and i](t,x) = w{x,t). 
Hence we conclude that the scalar super-field equation (ptt = 4>x + 4>t ■ T^4> 
provides a one-component representation for that system; we recall that a dif- 
ferent one-component representation of the fermionic extension for the Burgers 
equation derived in [6] is an evohition equation on an associative algebra. The 
geometry of the fermionic (supcr-)cxtcnsion (3) was extensively studied in [6]; 
in particular, recursion operators were constructed for its symmetries. We em- 
phasize that the correlation between system (2) and the extended Burgers 
equation was not noticed in [18], where (2) was first described. 

The geometry of system (2) itself is essentially nonlocal, which is illustrated 
by calculating the conserved currents using SsToOLS, trivializing them, and 
constructing the recursion operators that involve the new nonlocalities. It is 
very likely that system (2) has only one conserved current which is already in 
use for specifying 0. In the second layer, many nonlocal conservation laws and 
hence many new variables appear. This is discovered by SsToOLS as follows. 

The calculation of conservation laws for evolutionary super-systems with ho- 
mogeneous polynomial right-hand sides is performed by using the procedure 
ssconl: 

ssconl(N,tw,inincw,maxcw, afwlist .abwlist ,pdes) ; 

where 

N ... the number of superfields 

tw ... 2 times the weight 
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mincw . . . minimal weight of the conservation law; 

maxcw . . . maximal weight of the conservation law; 

af wlist . . . list of weights of the fermionic fields f (1) ,. . .,f (nf ) ; 

abwlist . . . list of weights of the bosonic fields b(l),. . .,b(nb); 

pdes . . . list of the equations for which a conservation law must be found. 



The ansatz for the differential polynomial components of a conserved current 
is composed in full generality (certainly, it has nothing to do with the axioms 
for the symmetry flows). Again, the global positive integer variable max_deg 
determines the highest power of a bosonic variable of nonpositive weight or a 
fermionic variable of negative weight and all their (super-) derivatives in any 
ansatz. The procedure ssconl (and wgts and linearize as well) is indifferent 
w.r.t. the presence of assignments = and => in pdes. 

The fact that the current is conserved on a given system pdes leads to the 
algebraic system for the undetermined coefficients, which is further solved 
automatically by Crack [16]. Having obtained a conserved current, we define 
the new bosonic or fermionic dependent variable (the nonlocality) using the 
standard rules (1). An example illustrating the run of ssconl is given in the 
test output. 

So, in addition to the potential 0, let us construct the fermionic variable v 
whose weight \v\ = | is minimal (other admissible nonlocal variables have 
greater weights): we calculate conservation laws involving and then we set 
Vt — Vh- (t)fh-\- fx(t)f and Vv — —Vh ■ fb+T>f •T>b-(j) + b^cf)/. Now there appear 
nontrivial solutions to the determining equations for recursion operators. By 
convention, the capital V denotes the 'linearized' counterpart of v and like- 
wise $ for (j). Recall that the correlation between the 'linearized' variables in 
the solutions TZ obtained by SsToOLS and the recursion operators R that act 
on the symmetries is as follows, see (6): the local variables F, B, and their 
derivatives denote the corresponding components of the fiows, and the lin- 
earized nonlocalities, e.g. $, act by the rules that follow from their derivation 
formulas; we thus have the interpretation $ = T>~^[B). 

We obtain the recursion of zero differential order with nonlocal coefficients: 




-Vb ■ (f)fB + (f)vF + v-B\ 
Vb(f)fF -vF + v(j)- B )' 



Also, we get another operator, which is nonlocal and has nonlocal coefficients, 

_ fVb V(f) - Vf VB (t)f - Vf Vb + VfVbF(l) + VfV + V(f)fb\ 
^[-2] - yT>BVb(f)f + VbV-VbF(f)fb + VfVbV(l)f + VfV(l)-Vfb J' 

it contains (p as well as $ and V. The coefficients of the recursions "^[-5/2] and 
7?.[_3] found for \sii\ — — 2| and \sr\ = —3 are also nonlocal. 



14 



3 A super-field representation for the Burgers equation 



Consider the system 

ft^Vb, bt^Vf + b'', (4) 

which has the unique set of weights |/| = |6| = |, \Dt\ — |, {D^l — 1. 
Equation (4) admits the infinite sequence (7) of higher symmetries fs = (l>^ , 
bs = (p^ at all (half-) integer weights \Ds\ > \. Also, there is another infinite 
sequence (8) of symmetries for (4) at all (half-) integer weights \Ds\ > \ of the 
odd 'time' s. 

System (4) is obviously related to the bosonic super-field Burgers equation 

b^ = bu-2bbt, b = b{x,t,e). (5) 

We emphasize that the role of the independent coordinates x and t is reversed 
w.r.t. the standard interpretation of t as the time and x as the spatial variable. 
The Cole-Hopf substitution b = —u~^Ut from the heat equation = Uu 
provides the solution for the bosonic component of (4) . 

Let us introduce the bosonic nonlocality w{x,t,6) of weight \w\ = by trivi- 
alizing the conserved form of the first equation in (4). Namely, we specify the 
derivatives of w by setting 

Vw = -/, Wt = -b. 

Note that the variable w is a potential for both fields / and b. The nonlocality 
w satisfies the potential Burgers equation = Wu + such that the formula 
w — Inu gives the solution; the relation / = —Vw determines the fermionic 
component in system (4). 

Now we extend the set of dependent variables /, 6, and w by the symmetry 
generators F, B, and W that satisfy the respective finearized relations. 

Ft = VB, Bt^VF + 2bB, VW = -F, Wt = -B. 

The linearization correspondence between the fields is / i— > F, 6 i— > S, and 
w I— > W. In this setting, we obtain the recursion of weight = —1: 



_ (F^-VfF + lJV\ 
"^^-'^-U-VfB + b^W ^ 



^D,-Vf + f,V-' ^ 
^ b,V-^ D,-Vf^ 



(6) 

The method for constructing 7?-[-i] is illustrated in the test output. 

The recursion 7?.[_i] provides two sequences of higher symmetries for sys- 
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tem (4): 



(ft] ( VK-VfVh-Uh \ //A (f.,~2Vff^ 

[bj ^ [Vf, - {Vff - b^Vf + bbj ^ ■ ■ ■ ' [bj ^ [b,, - 2VfbJ ^ ■ ■ 

(7) 

The same recursion proliferates two experimentally found first order flows to 
two infinite sequences of symmetries with the odd parameters s: 

(Vf\ fVU-{Vff-U\ / fVb-bVf + b, \ 

[vbj ^ [vb, -vfvb- bj) ^ ■ ■ ■ ' [bvb -fvf + f,- fb^J 

(8) 

Remcirk 4 System (4) is not a supersymmetric extension of (5); it is a 
representation of the bosonic super-field Burgers equation. The flows in (7) 
become purely bosonic in the coordinates Vf, b. The standard recursion 
R = Dt + ^b + ^bfD^^ for the Burgers equation acts 'across' the two sequences 
in (7) and maps {ft,bt) ^ {fx,bx)', again, we note that the independent coor- 
dinates are reversed in (5). 

However, from the above reasonings we gain two sequences of symmetries (8), 
which are not reduced to the bosonic {x, t) -independent symmetries of the 
Burgers equation. We finally recall that the Burgers equation (5) has infinitely 
many higher symmetries [19] that depend explicitly on the base coordinates 
X, t and thus exceed the set of axioms on p. 5. 



4 A multiplet of super-systems 

In this section we consider the systems 

ft = -afb, bt = Vf + b\ (9) 
Here the normal growth equation upon 

f{t, X- e) = /(O, X- e) ■ exp (-a |J 6(t, x- 6) dr) 

is coupled with the perturbed blow-up equation upon x; 9). System (9) has 
multiple homogeneities, and we let the tuple |/| = |6| = |, \Dt\ = |, \Dx\ = 1 
be the primary 'reference system.' The symmetry integrable cases differ by the 
values of the parameter a. We discover that for a — 1, 2, and 4 systems (9) 
admit the symmetries with both even and odd times s, s. The equations in this 
triplet demonstrate different properties. The geometry of the a = 2 system is 
quite extensive: this system admits a continuous sequence of symmetries for 
all (half-) integer weights \Ds\ > |, four local recursions (one is nilpotent), and 
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three local super-recursions. The equation for a = 1 admits fewer structures, 
and the case a = 4 for (9) is rather poor. All the three equations possess 
higher symmetries with parity reversing times s. In what follows, we indicate 
the presence of recursion operators for these systems in certain weights \sji\, 

Next, we find values of the parameter f3 = —1/a such that the respective 
systems (9) possess symmetries with parity reversing times s but have no 
higher commuting flows with the times s. Eleven cases are then realized and 
we specify the admissible weights |s| for each (3. 



4-1 Case a = 2. 



First we fix « = 2 and consider (9): we get ft = —^fb, bt = Vf + 6^. The 
weights for symmetries arc \Ds\ = |, \Ds\ = 1, and then equation (9) admits 
a continuous chain of symmetry fiows for all (half-)integer weights \Dg\ > 2^. 
Surprisingly, no nonlocalities are needed to construct the recursion opera- 
tors, although there are many conservation laws for this system. Hence we 
obtain purely local recursion operators TZ that propagate the symmetries: 
ifi = {F, B) ifi' = TZ for any ip. We get an operator of weight — 3|, 

^[-7/2] = {{pfffF + 6{Vfyfb^F + 12V f ■ fb^F + 8/66f) ' 

which is triangular since TZ'' does not contain B and, moreover, which is nilpo- 
tent. The above recursion is a recurrence relation [6] that is well-defined for all 
symmetries of (9). The recursion 7?.[_2] of weight — —2 is also triangular; 
we have 

^ ^VF Vff + 1 IV F fb' + I {VffF + 3Vf Fb^ + \UFf \ 

^[-2] = 

1 IVB fb^ + 8Vb Fb^ + 22Vb fBb + 7{Vf fB + 
\ UVf Bb' + fVf VBf + IV f VbF+'^b,Ff + f,Fb + 5UBJ 

Further, we obtain the recursion TZ[-5/2] of weight 2|; its components are 

1 = -2VbFfb'^ -VFVf fb-VF fb^ - y^Ffb-2Vf fBb'^, 



2tl " " " "2. 

-Rj; 1 = VB fb^ + VbFb^ + VbfBb^ + lVf^Ff + 

"^^^ + \VFb^ + \VF{Vff + VF Vf 6^ + IVF fj + {VffBb + 
+ Vf Bb^ + Vf VB fb + Vf Vb Fb + Vf Vb fB + |P/ F^f + 
+ \Vf f^F + \b,Ffb + \F,fb'' + If^Fb'' + If.fBb. 
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—3 is huge. 



For a = 2, system (9) admits at least three super-recursions ^{R-^, R'^) such 
that the parities of R'^ and R^ are opposite to the odd parity for / (and 
hence for F) and to the even parity of b and B. This is possible owing to 
the presence of the variable sr. The triangular zero-order super-recursions 
arc :^[_2] = ^'^fFfb + 8Ffb^, n\_^^ = -AVbFfb + 2{VffF + %Vf Fb^ + 
AVf fBb - f^Ff + AFb'^ + 8fBb^ and 

^ ^( -Vf^F-2^Fb' \ 

[-^^] [Vb UF - Vf b,F + Vf UB - 26,F62 + 2UBb^l 

for the weights \sr\ = —2 and \sr\ = — 2|, respectively; the third super- 
recursion found for \sr\ — — 2| is very large. Quite naturally, system (9) has 
infinitely many s-symmetries if a = 2. 



4.2 Case a ^1. 



Setting q; = 1 in (9), we obtain the system ft — —fb, bt — Vf+b"^. The default 

set of weights is again |/| = \b\ = |, \Dt\ = |, and \Dx\ = 1. The sequence of 
symmetries is not continuous and starts later than for the chain in the case 
a = 2. We find out that there are symmetry fiows if cither \Ds\ = \Dt\ = | 
(the equation itself), \Ds\ = \Dx\ = 1 (the translation along x), or \Ds\ > 3| 
such that a continuous chain starts for all (half-) integer weights \Ds\. 

Similarly to the previous case, no nonlocalities are needed to construct the 
recursions, which therefore are purely local. The recursion operator 'R-[_^/2] ~ 
0, 7^|'_5/2] = C^ffFf + iVfFfb'^ + \Ffb^ of maximal weight |s^| = -2\ is 
nilpotent: Ti? = 0. For the succeeding weight \sr\ = —3, we obtain a nilpotent 
local recursion 7?-[_3] whose components are given through 

7^f_3] = lVF{Vfff + IVFVffb' - liVffF - \{VffFb^ + 
+ bVf Vb Ffb + fVf f,Ff + ff.Ffb^, 

7^f_3] = Vf^ Ffb - ^VF Vb fb^ - ^-fVF VfVbf + IIVF fjb + 
+ \{Vf fVBf + l{VffVbF + \VfVBfb'^ + IVfVbFb'^ - 
- 55V f Vb fBb + f b,Ff + Vf fJB + fb.Ffb^ + '-ffJBb'. 

The differential order of 7^[-3] is positive. 
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4-3 Case a = 4. 



Finally, let a — A; then system (9) acquires the form ft — —^fh, bt — b'^ + Vf. 
Again, the primary set of weights is |/| = \b\ = |, \Dt\ = |, \Dx\ = 1. 
System (9) admits the symmetries {fs, bs) of weights \Ds\ = |, 1 or \Ds\ > 3| 
w.r.t. the primary set above. This situation coincides with the case a = 1. 
Again, no nonlocahties are needed for constructing the recursion 7^[_7/2] of 
weight = — 3|; this operator is relatively big: 

1 = -l2VbFfb'^ -VF{Vfffb-AVFVf fb^ -WF fb^ - 

'"^2l _ 4(pjyfBb'^ - AVf Vb Ffb^ - IVf f^Ffb - 12V f fBb^ - 2f^Ffb^, 
, = 3VBfb' + 3Vb Fb'' + Wb fBb" + \Vf., Vf Ff - \Vf,, Ffb^ + 

2 +VFVb fb^ + \VF {Vff + \VF {Vffb'^ + {VF Vf b^ + 
+ Vf fj + \VF fjb^ + {VffBb + 4{VfYBb' + 

+ {VffVbFb + {VffVbfB + liVffFJ + \{pfYf,F + 

+ 4r>/ VB fb^ + 4X)f Vb Fb^ + IQVf Vb fBb^ + Vf F^fb'^ + 

+ pf f.Fb-" + IVf f.fBb + 2b,Ffb' + FJb' + y,Fb' + f.fBb' + 

+ \VF b^ + VF Vf Vb fb + {VffVB fb + ?>Vf Bb^ + §P/ b^Ffb. 

No nilpotent recursion operators were found for system (9) if a = 4. 

Remcirk 5 We observe that an essential part of recursion operators for sti- 
persymmetric PDE are nilpotent. At present, it is not clear how the nilpotent 
recursion operators contribute to the integrability of supersymmetric systems 
and what invariants they describe or symptomize. We emphasize that this 
property does not always originate from the rule 'f-f — 0', but this is an 
immanent feature of the symmetry algebras. 

Problem 6 Construct an equation S that admits nilpotent differential re- 
cursion operators . . . | R^' =0} which generate an infinite sequence of 
symmetries </?, Rij^{(p), Ri^o R^^[(f), ... for S. Here we assume that at least two 
operators (without loss of generality, Ri and R2) do not commute and hence 
the flows never become zero. 

4-4 The systems with parity reversing parameters s 

Now we consider eleven systems (9) that admit symmetries with parity re- 
versing parameters s but do not possess any commuting s-flows except for the 

cases s = t, \Dt\ = ^, and s = x, \D^\ = 1. These systems are enumerated 
by the parameter f3 = —1/a whose usage makes the ordering simpler: we con- 
jecture that the systems possessing higher s-symmetries constitute an infinite 
family corresponding, in particular, to (half-)integer positive /?'s. In Table 2 
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we show the numbers of symmetries with a certain weight |s| for each of the 
systems (9) specified by (3. For convenience, we use the computer notation for 
the weights, that is, we multiply them by 2. Empty boxes correspond to no 
symmetries in that particular weight. 

Table 2 

The symmetry struct ure; for special s_ystcms (9). 



13 i 1 2\Ds\^ 
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6 
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5/2 














• 
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-3/4 
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-1/3 
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-1/8 


















• 


3 


















• 



5 A system with parity reversing times 

Let the time f and the parameters s be parity reversing. Then there is a unique 

system that satisfies the axioms and the weight assumptions |/| = \h\ — \Di\ — 
1 

ft = Vf + b\ k = Vb + fb. (11) 

Clearly, the weights are uniquely defined in (11). The numbers of symme- 
tries for system (11) are arranged in Table 3; note that the integer weights 
correspond to |s| and the half -integer values stand for 

Table 3 

The number of symmetries for system (11). 



\Ds\ e Z, \Ds\ e'L+\ 


1 

2 


1 


1^ 


2 


2i 


tJsym 


2 


4 


2 


1 


1 



The fiows commuting with (11) are proliferated by three local recursion oper- 
ators. The nilpotent super-recursion 7?.[_i/2] of weight = —\ and the local 
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recursion with \sji\ — —1 are given through 

V ( \ V (b-VB-Vb-B-Vf.F\ 

Another parity- reversing recursion for system (11) is 

^/ fVb-B + fVf-F-fb-VB 

y^bVb -B + ^bf-VF- . x>B + \bVf ■ F - fVf ■ B 

here we have \sr\ = — 1|- 

In the test output, we present the SsToOLS calls that demonstrate how the 
recursion 7^[-i/2] is obtained and we explain how the output is translated into 
the geometric language. 



6 Conclusion 

Using the presented computational techniques, we discovered that four classes 
of nonlinear coupled boson-fermion systems satisfy the set of axioms, which are 
believed natural in view of the nontriviality and nondegeneracy they ensure. 
Two of the classes are related to the Burgers equation, which is one of the 
most relevant and model-like equations of mathematical physics. Namely, they 
are a fermionic (super-) extension and a supersymmetric representation of the 
Burgers equation, respectively. The latter admits symmetries such that both 
the bosonic and fermionic fields are involved and such that the flows do not 
retract to the standard, purely commutative setting. Also, we foimd that a 
unique evolutionary system in this set of homogeneity weights has the parity 
reversing time t. Finally, we obtained a multiplet of coupled normal evolution 
and explosion systems; the multiplet seems to be essentially infinite containing 
all (half)integer values /3 G |M of the parameter at the least. Only three of 
the systems from this multiplet corresponding to a G {1,2,4} admit higher 
symmetry flows with the parity preserving parameters s; the parameters s 
for higher symmetries of other systems are parity reversing. We do not know 
whether these three systems are distinguished by any relation to the three 
classical sequences of the complex semi-simple Lie algebras. Hence we see 
that the suggested classiflcation scheme yields both generalizations of known 
integrable systems and also new models, which are expected to be relevant 
in nature in view of their vast symmetry properties: recursion operators were 
constructed for the above systems. 

The search for super-systems together with their symmetries and further re- 
vealing the recursion operators were performed by the SsToOLS package, 
which was used in two modes, respectively. The breadth flrst search method 
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for the recursions, see Remark 3, which is vahd for the systems that admit 
multiple homogeneities, proved to be a serious help in debugging SsToOLS. 
Thus we conclude that the general classification problem (see section 1.1), 
which is an immense task of modern computer physics, is becoming tractable 
with this computer algebra package. 
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TEST RUN OUTPUT 

Here we illustrate the process of computing symmetries, recursions, and con- 
served currents by SsToOLS. First we consider the weights |/| = |6| = \Dt \ = \ 
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and find all systems that satisfy these homogeneity weights and admit higher 
symmetries. We discover that system (11) is a unique solution obtained by the 
call 

ssym( 1 , 1 , sw , {1} , {1} , {} , {} , {} , {tpar , spar}) ; 

where sw G {1,2,3,4,5} and the flag spar is omitted whenever sw is even. 
The weights for this system arc uniquely defined, which is confirmed by the 
call 

wgts({df (f (l),t)=b(l)**2+d(l,f (D) , 

df (b(l),t)=f (l)*b(l)+d(l,b(l))},{l},{l},maxwt); 

here majcwt can be any positive integer. 

Now we present the intermediate computations that result in the recursion 
7^[_i/2] for (11), see (12). We compute the linearization of system (11), see 
section 1.2 for its definition: 

linearize({df (f (1) ,t)=b(l)**2+d(l,f (1) ) , 

df (b(l),t)=f (l)*b(l)+d(l,b(l))},l,l) ; 

The output contains system (11) and its linearization: 

df (f (2) ,t)=2*b(2)*b(l) + d(l,f(2)); 

df (b(2),t)=d(l,b(2)) + f(2)*b(l) + f(l)*b(2); 

The linearization correspondence between the fields is / = f (l) i-^ F = f (2) 
and b = b(l) B = b(2); the scaling weights of F = f (2) and B = b(2) are 
equal to the weights of / = f (1) and 6 = b(l) . Note that the numbers nf and 
nb of fields are automatically doubled. The hnearized system is calculated at 
most once in any run of SsToOLS; the output is pasted into the next call as 
is. 

Now is the time to find the recursion such that \sr\ = —\ and hence sw(= 
—2\sr\) — 1. The reason why the substitutions => are used in the original 
equations is explained in Remark 1 on p. 9: by construction, the recursions 
are symmetries of the linearized system, and the original equations are used 
for intermediate substitutions only. We thus have 

ssym(l,l,l,{l,l},{l,l}. 

{df (f (1) ,t)=>b(l)**2+d(l,f (D) , 
df (b(l),t)=>f (l)*b(l)+d(l,b(l)), 
df (f (2) ,t)=2*b(2)*b(l) + d(l,f(2)), 
df (b(2),t)=d(l,b(2)) + f(2)*b(l) + f(l)*b(2)}, 

{},{}, {tpar , spar , lin}) ; 
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The output contains a unique solution: 



1 solution was found, 
df (f (2),s)=f (2)*f(l); 

df(b(2),s)= - l/2*f (2)*b(l) + f(l)*b(2); 

This is precisely the first recursion in (12) for system (11). The result is 
achieved by the intermediate call of Crack [16] from ssym; the solver can 
be handled in full automatic mode by typing g 1000. 

Let us present a similar calculation that yields the local recursion (10) for 
system (9) with a = 2; note that this time the weights are not uniquely 
defined and we thus have the multiple homogeneity case, see Remark 3. The 
input is 

hom.wei : ={{16 , {3 , 3} , {2 , 2}}} ; 

ssym(l,l,7,{l,l},{l,l}, 

{df (f (l),t)=>-2*f(l)*b(l), 

df (b(l) ,t)=>d(l,f (l))+b(l)**2, 

df (f (2) ,t)=-2*f (2)*b(l)-2*f (l)*b(2) , 

df (b(2),t)=d(l,f (2))+2*b(l)*b(2)}, 
{},{},{lin, filter}); 

the output is the recursion (10): 

1 solution was found, 
df (f (2),s)=0; 

df (b(2) ,s)=d(l,f (l))**3*f (2)*f (1) 
+ 6*d(l,f (l))**2*f (2)*f (l)*b(l)**2 

+ 12*d(l,f (l))*f (2)*f (l)*b(l)**4 + 8*f (2)*f (l)*b(l)**6; 

Note that the admissible set of weights |6| = \B\ — 0, |/| = \F\ = — | for 
system (9) and the value = 1 would require the upper bound max_deg : =6 ; 
for this solution. 

By shifting the primary weight sw and/or the additional values sw_i in hom_wei, 
one can inspect how new recursion operators appear and vanish. 

A computation of a conserved current using ssconl is done as follows; here 
we get the conserved current for system (2) having already extended it by the 
variable = f (2) , see p. 13. The input is 

ssconl(l,l,5,5,{l,0},{l}, 

{df (f (l),t)=d(l,b(l))+f (l)*b(l), 
df (b(l),t)=d(l,f (D), 
df (f (2),t)=>f (1), 
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d(l,f (2))=>b(l)}); 



Then the output produced by SsToOLS is 

NEXT: BOSONIC CONSERVATION LAWS OF WEIGHT 5 

>»» Non-trivial conservation law: 

Pt = Db*f (l)*b - Df (l)*Db*f(2) - b *f(2)*f (1) 

X 

Pd(l) = Db*f (2)*f (l)*b + f(l) *f(2)*f(l) 

X 

By trivializing the conservation relation Dt(Pt) +'D(Pd(l)) = 0, we construct 
the fermionic variable v of weight l^;! = |, see section 2. 
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